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CFD  METHODS  FOR  DRAG  PREDICTION  AND  ANALYSIS  CURRENTLY 

IN  USE  IN  UK 

by 

P.  R.  Ashill 


SUMMARY 


Computational  methods  developed  in  UK  for  the  prediction  of  the  drag  of 
aircraft  components  at  subsonic  and  supersonic  speeds  are  critically  reviewed. 
In  many  cases,  the  flow  modelling  isfound  to  be  lacking  in  certain  respects. 
Despite  this,  however,  the  review  suggests  that  these  methods  have  a  useful 
function  both  in  the  early  stages  of  aircraft  design,  when  they  may  be  used  to 
study  differences  in  the  drag  of  various  shapes,  and  later  in  support  of  wind- 
tunnel  tests  as  a  diagnostic  tool  and  also  to  'extrapolate'  the  data  to  'full 
scale'.  i. 
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LIST  or  SYMBOLS 


wing  ratio 

local  atraamwla*  chord 

<1rag  co<»rriclcnt  based  on  wing 
planform  area 

local  drag  roafricletit  based  on 
local  ctiord 

drag  ro»rriciant  baaed  on  aurface 
area 

cowl  preaaitre  drag  coefficient 

(PU  ?") 

notional  drag  coefficient  per 
aiirftce  • 

drag  coefficient  based  on  frontal 
area 

akin  friction  coefficient  baaed 
on  free  stream  dynamic  preasure 

lift  coefficient  based  on  wing 
planform  area 

pitching  moment  coefficient  baaed 
on  wing  planform  area  and  mean 
chord,  nose  up  positive 

preaaiire  coefficient 


free  stream  Mach  number 

Mach  number  of  the  flow  component 
normal  to  and  Just  upstream  of  the 
shock 


Reynolds  number  baaed  on  atreamwlae 
chord 

wing  planform  area 

plane  normal  to  free-atrtsui  vector 
and  downstream  of  aircraft 

dlatanee  normal  to  the  wing  surface 

angle  of  Incidence 

incremental  part  of 

non-dlmenalonaX  apanwiee  dlatanee 


balance  mtaaured 
body  alone 

akln-frletlon  component 

normal-praaaura  component 

tralllng-vortei  component 

vlecoua  or  boundary-layer  oomponent 

wave  component 

notional  wing  alone 

far  upstream 


I  IMTRODUCTIOM 

In  the  Technical  Evaluation  Report  of  the  ROARD  Conference  on  'Aerodynamic  Drag' 
held  at  Icmir,  INirkey  in  April  1973*  it  wae  concluded  that  *a  aommrthenalft  drag  predic¬ 
tion  methodp  valid  for  the  main  claaaea  of  aircraft  and  baaed  entirely  on  theory,  la  not 
likely  to  be  poaalble  for  a  long  time  to  eeae**  fifteen  yaara  later,  the  wholly  theoreti¬ 
cal  prediction  of  aircraft  drag  to  a  aatlafactory  atandard  of  accuracy  la  atill  not 
poaalble.  However,  this  period  haa  seen  eomiCerabla  progrtae  in  tha  davalopaant  of  flow 
algorithms,  notably  for  trenaonle  flows,  and  a  roduatlon  In  tha  eoat  of  eomputationa  of  at 
least  two  orders  of  isagnltude’. 
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These  developments  hsve  encoureged  the  Increasing  use  of  CFD  In  the  design  of 
eircreft  from  the  preliminary  stages,  through  the  deeelopnent  phase,  to  pre-production. 

In  the  early  atagee,  approximate  CPD  methods  (eg  Ineiaeid  methods)  provide  the  project 
engineer  with  simple  tools  for  selecting  suitable  designs.  Later,  during  the  development 
phase,  increased  reliance  is  placed  on  more  complex  CFD  methods,  including,  for  example, 
viscous  effects.  Combined  vith  data  from  csrefull/-eondueted  eind-tunnel  testa,  these 
methods  enable  the  designer  to  diagnose  sources  of  excess  drag  and  to  predict  the  drag  of 
modified  shapes,  tfsed  In  this  «ay,  the  methods  need  only  be  reliable  in  their  predictions 
of  small  drag  dirferencea  and  thus  it  is  not  necessary  for  the  floe  modelling  to  be  pre¬ 
cise  so  long  ns  the  main  featurea  of  the  floe  are  repreeented.  At  thie  stage  CfD  alto  has 
an  Important  supporting  role  in  the  elnd-tunnel  testa  for 

(i)  Eatabllahlng  a  basis  for  simulating  full-scale  flows  In  the  eind  tunnel  and,  ehere 
necessary,  extrapolating  the  tunnel  data  to  full  aeale; 

(11)  Calculating  tunnel  eall  and  model  support  Interference. 

Although  the  second  application  la  Important  It  la  indirect  and  is  not  considered 
further  In  this  paper. 

Finally,  before  production,  it  la  neeeaaary  to  guarantee  perfemanee  predictions 
from  prototype  night-teat  data,  and.  In  this  phaae,  CFD  has  a  poaaible  role  in  the 
interpretation  of  th#  flight  teat  data.  Again,  however,  thie  septet  la  not  dlaeuaaed  in 
the  paper. 

Thie  paper  reviews  current  UK  CFD  methods  for  drag  prediction.  Where  possible,  the 
predictions  are  compared  with  measurement^  otherwise  results  of  csleulstlona  are  included 
to  llluetrate  the  use  of  the  methods  In  aircraft  design.  Because  of  llmltstlona  on  the 
length  of  the  paper  the  review  Is  not  exhaustive  but  it  is  hopad  that  the  paper  gives  the 
flavour  of  UR  activities  in  this  field. 

Following  a  discussion  of  general  aspects  of  drag  prediction  in  section  2,  the 
paper  reviews  mwtltods  for  subsonic  aircraft  In  section  3  and  for  supersonic  slrersft  In 
section  A. 
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OKNKKAL  CONf^IllFHATIONS 


Two  alternative  procedures  are  available  for  obtaining  drag  from  CFD  predictions, 
as  shown  in  Fig  I;  t»ie  first  or  'local*  method  involves  integration  of  the  atreamwise 
contributions  of  the  forces  due  to  normal  pressure  and  akin  friction;  the  second  la  a 
'field'  method  re<julrlng  an  Integration  over  a 
plane  normal  to  th*  free  stream  and  downatream  of 


the  aircraft. 


•T' 


HctM  1  toCAt 


to 


The  susceptibility  of  the  'local*  method 
truncation  errors  Is  well  known  and  results 
obtelned  by  this  technique  should  always  be 
checked  for  the  effect  of  grid  spacing.  The 
'field*  fsethod  may  alao  be  aenaltlve  to  grid  den¬ 
sity  but,  aa  yet,  there  la  little  experience  on 
which  to  base  a  Judgement  of  this  procedure. 

Inveatlgatlng  the  drag  of  an  aerofoil 
inferred  from  calculations  by  an  Inviscid  Euler 
code,  Tu  et  el'  showed  that  both  the  'local*  and 
'field*  methods  Incorrectly  gave  non-*ero  drag  for 
a  eubcritlca]  flow.  Lock^  attributed  thla  problem 
to  the  generation  of  spurious  entropy  near  the 
leading  edge.  Tiius  It  would  appear  that  further 
development  of  flow  algorithms  la  needed  before 
the  'field*  method  can  be  used  with  confidence. 

Chi  the  other  hand,  with  possible  enhancements  In 
mind.  It  may  be  noted  that  the  'field*  method, 
unlike  the  'local*  method,  does  not  depend 
directly  on  details  of  the  aircraft  geometry  and 
My  thus  find  an  application  to  the  prediction  of 
the  drag  of  complex  configurations. 
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With  the  pl.ne  ’T'  t*ken  «urriel*ntl7  f.r 
down. treat.,  the  earloue  tenae  In  tha  'ftald' 

Integral  tmt  be  e.panded  In  poverc  of  the  pertur¬ 
bation  eeloeitiaa  (non-dinienaionalia|d  uith 
reapeet  to  frae-atraan  apeed).  I,oak  ahouad  that, 
to  an  order  of  approii.ation  that  it  adequata  for 

aubaonle  tranaport  aircraft  at  eruitt  eonditiona,  thit  aaprataion  reducaa  to  tha  elaaaioal 
'far  fi.ld'  integr.l  which  can  be  divided  into  three  eoegtenente  aa  ahonn  in  Pig  t. 


ri,.1  Iwa  Mtkedi  et  dcIcrMnki,  dray 
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Lock*  observed  thst  the 
dreg  co«sponenes  of  vjitge  eoaid 
be  determined  most  eonvenlentljr 
end  eeeurstelj  by  relstlng  flow 
conditions  st  *T*  to  those  on  or 
nesr  the  wing.  The  three  drag 
eosiponento  are  treated  as  rolloea: 

(a)  Have 


ohac 


voaux 


VISCOUS 


HAVC 


On  the  reasonable  assunp^ 
tion  that  the  flow  downstream  of 
all  the  shocks  Is  Isentroplc  and 
adiabatic,  wave  drag  la  determined 
bj  the  reduction  In  total  pressure 
across  each  element  of  the  shock 
ayaten.  This  statement  has  no 
meaning  for  potential  flows  but 
methods  have  been  developed  In  UK 
for  Inferring  wave  drag  from 
potentlai*riow  solutions.  A 
method  for  aerofoils  at  subsonic 
free<'8tream  speeds  due  to  Billing 
and  Bocci),  which  has  led  to  the  development  of  the  computer  program  known  as  RACRCONT, 
relates  each  element  of  the  shock  to  s  Bsnklne-flugonlot  shock  of  the  same  atrsngth,  le 
having  the  same  Mach  number  normal  to  and  Just  upstream  of  the  shock,  H||  .  Billing  and 
Bocel  also  asaumed  that  the  local  flow  la  normal  to  the  shock.  This  assumption  la  reason¬ 
able  for  inviscld  Dows  at  high  subsonic  speeds  but.  In  vlseous  flows,  where  the  Inter¬ 
action  beween  the  shock  and  the  boundary  layer  causes  the  shook  to  be  oblique  nesr  the 
aerofoil  aiirface,  th«  method  probebly  overestlmatea  wave  drag. 


flf.2  Far  lleW  ssslyth  el  drsp  sf  wkit<  si  tsbtsak  ipe*d« 


This  approach  has  been  generalised  to  wing  flows  by  AllwrightS  except  thst.  In  his 
method,  no  saaumptlona  are  made  about  the  direction  of  the  flow  Just  upstream  of  the 
shock . 


In  ose^e  wher^  details  of  the  flow  field  are  not  known  or  a  rapid  Indication  of 
wave  drag  Is  needed,  a  simple  method  due  to  Locks  Is  useful.  In  Its  two  dimensional  form, 
Lock*8  approach  is  similar  to  thst  of  Billing  and  Boecl  eaespt  that  It  uses  the  assumption 
thst  the  shock  wave  lies  along  the  normal  to  the  aerofoil  section  contour.  With  this 
assumption  and  by  retaining  only  the  first  term  In  the  Naclaurln  expansion  with  respect  to 
distance  from  the  serofoll  contour  for  the  gradient  of  ahoek-upatraam  Haeh  number  Hn 
normal  to  th^  aerofnii  contour,  Lock  obtained  the  following  tkpreaalon  for  wave  drag 

^  -  co,e .  Ojlia  /i_iju29!Y  -  »*(?  -  , 

*•  \  *'  /  |,,oU  ♦  0.?H,o) 

H«r«  fl  1.  fre.-str...  Mach  nunber,  k,  1.  the  loo.l  euretture  or  the  eeroroll 
section  et  the  foot  or  the  ehoek,  derlned  b,  the  .orrii  0  ,  snd  ,  1*  rree-*tr.M  djn.- 
iilc  preeeure. 

Equation  (U  Inpllee  that,  for  a  (lean  ealiie  of  Nmo  •  .action  raee  dra,  ]n  Lock', 
•pprexlnatlnn  depend,  onljr  on  the  local  radtua  or  eureatura  l/k«  •  Thla  la  an 
appropriate  length  ecale  an  long  aa  either  (a)  the  aeroroll  eureature  ehangea  aloe]/ 
upatreai.  or  th.  .hock  or  (b)  the  height  or  the  aboek  penatratlon  Into  the  Field  la  anall 
eotapared  ulth  .  Thua  for  Hlnga  elth  both  a  aurraca  eureature  that  ehangea  rapldljr 

Kith  acreamelee  distance  and  a  strong  shock.  Lock's  Mthod  M,  be  expected  to  glee 
Ineceurst.  predictions  or  eavs  drag  (ass  aaetlon  3.2). 

Sines  l.aek's  aiethod  utllltsa  ths  assuaiptton  that  ths  shock  la  norasl  to  the  asro- 
roll  contour  and  Is  based  on  «lng  surracs  eureature.  It  does  not  Include  the  erreet  or  the 
elecous/lnvlaeld  Interaction  betsssn  ths  boundary  Isvar  and  tha  shook. 

Lock  nodirsed  equation  (1)  to  allou  ror  ulng  suesp  by  using  tha  assuaptlon  that,  at 
each  «lng  section,  the  riou  la  Identical  to  that  oeer  an  Inrinlta  yaved  ulng  haelng  tha 
aasie  aneep  aa  the  shock. 

The  dctenalnatlon  or  «aee  drag  rren  solnttens  of  tha  tular  aquations  la  lass 
stralghtroreard  than  It  Tlrat  appaart.  At  notad  aboea,  apurloua  antropy  la  Inearltbly 
produced  upatressi  or  the  shock  rroai  areta  such  at  the  ulng  leading  adga  ahara  thtra  are 
rapid  changes  In  shape  along  tha  alng  chord.  Thua  aaea  drag  ealeulatlena  bated  on* the 
field  aiethod  can  be  slgniricantly  In  error.  Attaavta  to  infer  the  oaee  drag  tnm  tha 
antropy  rlae  teroos  the  ahock  are  eoeiplleated  by  nuaorieal  arrora  In  the  roglon  of  the 
ahock.  Methods  or  dealing  with  this  problea  haec  been  dlaoutaed  by  Sallal  and  LoekA, 

(b)  Vortex 

In  order  to  base  any  retaonabla  proapeet  or  ealoulatlng  this  eonponent  dlractly,  It 
la  neeeaaary  to  Ignore  the  rolling  up  or  the  tralling-eortea  sheet.  Considerable 
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•iMpliricatlon  iB  also  possible  if  the  downward  Inclination  of  the  ahaet  la  Ignored,  the 
resulting  expression  being  the  classical  contour  Integral  around  the  vortex  trace  In  the 
Trefftt  plane.  This  approach  la  probably  adequate  for  high  aapect-ratlo  elngs  at  low  to 
Moderate  llft^  <  0.5)  but  for  low  aspeet-ratlo  einga  at  hl^  lift  It  auat  be  of 
questionable  accuracy. 

(c)  Viscous 

In  two-dimensional  flowa,  viaeoua  drag  May  be  Inferred  froM  the  solution  for  the 
viaeoua  wake  far  downstream  but  this  would  not  aeea  possible  for  flows  over  finite  wings 
because  of  eompllcattona  arising  from  wake-edge  conditions*.  Therefore,  for  wings,  or  if 
an  accurate  solution  la  not  available  for  the  viscous  wake  In  two-dimensional  flew,  an 
extended  version  of  the  Squlre/Toung  formula  allowing  for  compressibility  and  wing 
sweeps**  may  be  »is*d. 

Unless  otherwise  stated,  the  *rsr-fleld*  method  la  used  in  drag  predictions 
discussed  later.  As  shown  in  section  3.3,  this  simple  framework  for  analysis  appears  to 
be  lustlfied  for  subsonic  transport  aircraft  at  cruise  conditlMis.  For  flows  with  power¬ 
ful  interactions  between  the  viscous  shear  layers,  the  shock  waves  and  the  trailing  vor¬ 
tices,  a  decomposition  of  this  kind  is  no  longer  valid  and  the  scope  for  diagnostic 
studies  accordingly  limited,  furthermore,  overall  drag  would  then  have  to  be  calculated 
using  either  the  'local*  or  field*  methods  with  all  tM  difficulties  that  Implies. 

3  METHODS  POn  SUBSONIC  AIRCRAFT 

3.1  Aerofoils 

Methods  for  aerofoils  are  viewed  in  UK  as  a  first  step  towards  the  development  of 
aatlsfaetory  flow  algorithms  for  wings  and,  as  such,  have  been  used  to  test  Ideas  on 
various  aspects  of  flow  modelling.  However,  aerofoil  methods  have  progressed  to  the  point 
of  being  powerful  denign  tools  In  their  own  right  and  are  currently  used  for  tasks  such 
ss: 


il  scHmS ,, 
pin  kvUV 
•r  Crsfi'w 
osfM  ft 
•erfmi  vskti 

U  setlmd 


(I)  selection  of  wing  sections; 

(II)  design  of  flaps  and  slats;  and 

(ill)  extrapolation  of  tunnel  data  to  'full  aeale*. 

The  majority  of  the  methoda  currently  in  uae  in  UK  (Fig  3)  are  of  the 
viacoue/lnviscld  Interaction  type  in  which  calculation  of  the  two  parte  of  the  flow  is 
performed  Interactively  end  Iteratively  to  numerical  eonverganee.  A  number  of  numerical 
schemes  are  used,  namely  Direct  (which  is  only  suitable  for  attached  flew), 

Semi-Inverse  (SI)  ^whjch  may  be  used  for  separated  Lsw-tateS  atUmSi 

flows)  and  Ouasl-Slmultaneous  (QS)  (which  is  AICORIThns 

equally  effective  for  both  separated  and  attached  C(K)C  ORCMAIOR^  WVISCD  VBCDUS 

flows).  Full  details  of  these  schemes  are  given  '  ‘J ~  _ _  m  . _ _ _ 

In  the  review  by  lock  and  WllUamsV  Sssrct  Matl  -IM—  ll  uctim 

In  the  remainder  of  section  3,1,  the  *****  ****  — MkKi)— 

methods  summarised  In  Fig  3  are  reviewed,  methoda  Vflmsiw  srcrsti' 

for  low  speed  (and  high  lift)  being  coneidered  In  ssfM  t 

section  3*  1.1  and  techniques  for  high  subsonic  Hliu  %Nt  •mwtwi  tcF*  wsrimi  < 

•p,eds  In  .ectlon  M.?.  U  Mtta 

3.1.1  Low  .D.ed 

UK  method.  Tor  e.lcuKting  drag  .nd  Milmum 
lift  or  ..rofolla  St  low  rre«-gtr..m  ap««d.  my  b«  u  Nnii.tMid  ■tUwdi 
•  ummrl.ed  ta  follnwa:  - -* —  - 

visiRui  Mam  g  fm. 

>7  UHW  Sf  IrM.  H  SSlW 

(1)  SIVP  (Swl-lnv.r.t  VUeou.  Program)  — "  mw  ml 

mlmatim.  .... 

This  method"’  la  reetrleted  to  alngle  aaro- 
folia,  and,  aa  Its  nama  auggeata,  utlllaaa  an  SI 

aeheme,  with  a  aurface-alngularlty  taehnlqut  for  _ 

the  Inelecld  flow  and  Integral  methoda  for  the  It  meim 

ahear  layer. .  The  turbulent  boundary-layera  ara  t**  mhH  el 'S"** 

calculated  by  the  lag-entralnment  mthodO  while  nHiinil 

the  laminar  layers  arc  computed  using  a  .n  . 

eompreaalble  eerslon  of  Thwaltts*  mathodl*. 

Further  allowance  la  mde  In  the  turbulent  ama.  Wn.  immnti 

boundary-layer  calculation  for  the  afreet  on  akin  g  wmti”  —W—  Ua.li 

friction  of  low  Reynolds-numbar  (It  a  local  ealua  ■ 

of  momentum-thlekn.aa  Reynold,  numb.r,  R.  ,  1...  ^  yt.  ...m 

than  about  1000).  However,  R,  la  not  alfowad  to  *2; 

fall  below  3?o  Ju.t  downatream  of  transition,  ,  . 

since  this  a  natural  limit  for  a  fully-dcealopad  -gSb^ 

turbulent  boundary-layer.  In  addition,  the  seeon-  ll  •  li.  I.himmi.i 

dary  Influence  of  flow  curvature  on  turbulence 

structure  1.  Included  In  the  'lag'  equation.  Pi,.]  UK.  CPD  mtlhadi  far  .tr.l.a  Ra,  prcdkllan 


cooc 

ORIGMAIORI 

SWF 

vamst* 

.  NlOA 

Ntvmif  * 

Ntftr  1 

Wssiu 

FClIM 

SS.:.. 

visiRUil  Pkwm  a  I  fam.  rntlirmt  ,i 

I  tmn  II  »Fi.  »>  'rm>-  — — 


CFher  a  a.iii  imimlwt  U  ■ 

im  !•  Ml  mPfirtiwi.mMcii— 

hR.  1  M 
pmaiiii  imi. 


a  u.H«"  •'*  — 

om.  p.m»  **— 

a  arm  »  nil..  t(  llltr  _ IM  — 


7 


flnuMy,  Btan'inr'l  «har»-param»t.er  ralatlonahlp*  *  !•  replaced  by  one  that  la  aiore 
aultuMe  for  aeparet^d  flowa*.  No  allowance  is  made  for  the  'higher-order*  effects  In  the 
alr*am«Ioe  woiawntum  eTuatton  due  to  normal  pressure  gradients  and  Reynolds  normal 
{^trexars.  The  latter  *hlgher-order'  effect,  which  la  the  more  important  of  the  two,  la 
not  Included  because  correlationa  of  It**  are  of  doubtful  validity  for  flows  with  eiten- 
nlv*  r^glona  of  reparation. 


Whll*  the  method  giver 
close  predictions  of  sMxlmum  lift, 
as  Illustrated  for  two  different 
aerofoils  In  r*K  H,  ft  predicts 
much  lower  drag  than  lhat  sieasured 
on  the  aerofoil  na(W)-?  (Fig  5). 
This  discrepancy  mlaht  h» 
eiplatned  by  results  of  calcula¬ 
tions  which  sugr,*st  that  the  tran¬ 
sition  trips  us«d  1»>  this 
esperlment  wer»  not  adequate  over 
the  entire  range  of  incidences 
tested*.  The  neglect  of  the 
Reynolds  normal-stress  t*rm  aien- 
tloned  sbov*  may  a  loo  h» 

s 1 pnl r 1 csnt . 
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The  viscous  'pacVage'  in 
this  pr'^Ki  am  has  been  written  so  Fig  4 

thsf  It  can  readily  b«  coupled 
with  other  1nvlar1<|  methods,  and 
It  haa  al"o  been  used  In  tUe 
FFl.HA  and  Rrltlsh  Aerospace  (PAel 
Filler  codes  des«'r1hed  lat^r. 

fill  (Klgli  Lift  Design  and  Analysis) 

pwvoir>p«/i  f  .>  #*alculate  flows  over  multi-element 
sor'«f‘'>i  I  r ,  ihin  melhod''*  uses  the  Hlrect  coupling 
of  *  ho  oirMet*  MAV!.''«i*  (Multiple  Aerofoil 
VI  scour.  Iterative  ."ystom)  program  hut  haa  an  Improved 
su I  I  ng'i  I  a  I  1 1  y  motho«l  for  the  ( I  ncomppess  1  nle) 

1nvl-«'id  Clow’*.  Ar  ]n  r?ivr,  tlie  turbulent  boundary- 
Iryop-  and  1so|nt«d  wak^s  are  calculated  by  the  lag- 

•  nlralnmcnt  molhod.  l|o  allowance  Is  made  for 

*  h  1  fhof-oi  do  r  •  ^froofs  In  ilie  streamwla*  momentum 
equation  bu*  c  cn<  f.^ot.  io(i  fer  the  Influettce  Of  low 
n»yno]dc  nnmi'or  /.u  iur»'ui«nt  sk  1  n-f rl ct lon  Is  included. 


Vtrishea  of  Nfl  (ocfficlcaf.  C^.  with  ssQlt  of  kititftact.  o. 
for  two  icrefoHs 


FigS  Vorlillen  ol  Co  vilh  Ct  for  CAIW)*? 

wf.i.ii,  n  t  V.3  > 


♦  he 


Merging  (»i  the  wakes  from  ui'stream  elements  with  boundary 
iniogpsi  mei  iv'd  f'f  Irwin'’  snd  more  recently  by  a  method  due 


layers  Is  calculated  by 
to  Cross ’ ■ . 


51noo  (he  Direct  scheme  Is  used,  the  method  falls  where  separation  occurs  and  thus 
buid'io  sop.r i*at  I ons  occurring  In  re-entrant  or  *cove*  regions  are  empirically  modelled. 


rreflloMnnc  of  lift  and  drag  for  a  three- 
elcmeiii  sol•o^f^1  (  aio  shown  In  Fig  u.  Ihe  vlscous- 
Indu"e»f  In  MCt  fs  well  predloted  fop  angles 

of  Jin'idonco,  ,  up  to  ?n*  hul  ,  at  higher 
angles,  tiio  (|nw  separates  on  the  main  aerofoil 
and  consequently  the  method  fslla.  In  Ref  9  it  is 
argued  that  the  goori  agreement  between  calcula- 
tlofis  and  measurement  at  s  -  20’  la  to  some  estent 
fortiiltoiie,  the  lift  on  the  main  aerofoil  being 
overestimated  while  the  lift  on  the  other  two  ele¬ 
ments  Is  underestimated. 

The  estimates  of  drag  are  far  less  sstls- 
farfory  eaperlslly  as  the  stall  Is  approached.  Aa 
well  as  the  omission  of  'higher  order*  effects 
referred  to  above,  possible  reasons  include  the 
laek  of  pempreas 1  hi  11 ty  effects  In  the  calculation 
of  the  Invlaeld  flow  and  the  Inadequacy  of  the 
modelling  of  the  aerofoil  wake  In  the  region  of 
high  f low-curvatur*  above  the  flap. 

(Ill)  FELMA_y^nlt e  Element  Multiple  Aerofoil) 

Aa  Implied  above,  compressibility  can  eaert 
a  algnlflrant  influenoe  on  low  speed  flows  over 
multiple-element  aerofoils  at  high  lift  psr- 
tlriilarly  wliere  the  flow  accelerites  to  high 
speeds  locally,  eg  at  the  leading-edge  slat. 
FEI.MA'*  rwpi-osenls  comppean Iht 1 1  ty  In  the  tnviacld 
Daw  hy  solving  tho  eiact  potential  equation 


«  ?0* 


—  Inviscid 
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o  (lyrriment 


Fifl  Vsristfos  sf  lifl  with  snglr  of  intidrmt  ind 
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num»pic^lJy  by  a  rinl  technique.  As  noted  prevlousljr.  the  viscous  shear  layers 

ar<*  by  Utr  method  used  In  ??ivr  but.  In  contrast  to  HILDA  or  MAVIS,  FELMA  does 

not  represent  the  merging:  between  wakes  and  boundary  layers.  The  option  Is  provided  to 
use  either  SI  or  O^S  crtupllnae,  allowing  flows  with  separation  to  be  calculated.  Of  the 
two  schemas,  QS  Is  the  mor*  efrielent,  being  faster  than  Si  and  not  requiring  a  switch  from 
Direct  coupllnK  f-^r  the  ettac.hed  portion  of  the  flow  to  SI  coupling  In  regions  of 
separation. 

Cowparlaona  of  predictions  by  FELMA  snd 
■eaatirem^nta  of  lift  sort  drag  are  shown  In  Fig  7 
for  tha  MLh  t\0\  a** rofoi  1 /risp  configurations  1 
and  ?,  having,  r*ap«rtlvaly,  flap  gaps  of  2.6% 
and  1.3%  basic  aei-AroJl  chord.  The  calculation  of 
(sailmnin  lift  la  In  r«aaonaMe  accord  with  measure* 
went  for  tha  lergrr  of  th*  two  flap  gaps  but  for 
the  smaller  gap  th*  maximum  llTl  la  overestimated, 
possibly  beransp  an  observed  Interaction  between 
the  a»roroll  wskp  an<l  the  riep  boundary  layer  Is 
not  rapreaant^d  In  FELMA, 

Whlla  some  encouragement  can  be  drawn  from 
the  drag  pred\rt  in  Fig  7,  U  should  be  noted 
that  the  Hi,n  cnnT  I  gnrat  1  nna  are  somewhat  Idealised 
In  tha»  they  •]r>  n-'t  repreacriL  a  'cove*  on  the  main 
aeremu.  It  remalna  to  be  aeen  if  PEl.MA  offera 
improved  aeruracy  r^vci-  ti>at  of  HILDA  for  more 
prartlcal  cfitirig'uat  ton?!  where  the  merging  of 
wak*s  from  upatre»m  elements  and  boundary  layers 
may  be  art  Important  feature  of  the  flow. 

Overall,  the  ptesent  altnatlori  in  UK  as 
regarda  the  p 1 1 1. 1  on  of  drag  of  hlgti-llft  Sero- 
fejls  I »!  not  a  1  ♦■'’Pici  her  satisfactory.  TTtepe  gpe 
reasons  t ''  beHevo  ihat  this  arises  because  of 

In  the  mo.ioinnp  of  the  wake  of  the  main 
aeier^jl  fn  the  leptrvn  ^f  hlgti  f  1  ow-eti  rv8  ture 
above  the  Map.  }n  this  reglott  both  streamwlse 
ar>d  erosswjse  prernn-e  gradients  are  large  and 
hence  th«i  flow  there  Is  highly  ejitptlc  in 
cbsraei.ef.  TI»ms,  In  order  te  achieve  the  required 
aez-nracy,  it  may  l-e  necessary  to  use  one  of  the 
new  gencratloti  or  metl»ods  for  solving  the 
peynn  1  ds-aver*e.i ,  hs  v  1  •  PS  -  *  t  Ok  e  s  equation^*'.  How¬ 
ever,  these  meti'C'is  will  only  be  able  to  provide 
the  necessary  accinscy  If  turbulence  models  are 

wIJ'i,  eie  “  »' 1  gh  1  y -c  •  P  ved  wst^ssl'. 

lil£h_«ieed 

Pecause  of  the  Importance  of  being  able  to  estimate  accurately  aectlon  drag  for 
transpert-type  wings,  emphasis  has  been  placed  In  UK  on  the  development  and  validation  of 
t  ranscni  c- f  1  ow  co.i#»*s  (Pjg  1),  M-theds  currently  favoured  Include  those  based  on  the 
assumption  that  iiic  ittvlsrld  flow  is  potential  and  others  In  which  the  Euler  equations  are 
used  ♦c  simulate'  the  *oijtep  *  flow. 

using  potentlal-f iow  approximation 

The  code  V'^K’’  has  been  the  mainstay  of  wing  section  dealgn  and  analyala  In  UK  for 
over  years,  having  superseded  the  transonic  small'perturbatlon  code  VISTHAN^*.  VOK 
couples,  itt  the  iMrect  way,  a  numerical  solution  of  the  full-potential  equation  with 
Intagral  mcfiiodc  for  the  shear  layers,  the  laminar  and  turbulent  layers  being  calculated, 
respectively,  by  'niwaltefl  method”,  extended  to  allow  for  cosipreaalbillty,  and  the  lag- 

ent  ra  1  nmen  i  me  f  ho.i  •  '  . 

In  general,  VOK  gives  eatlafactory  predlctfona  of  drag  for  attached  flowa  but, 
wlicre  flow  ncpBration  Is  approached,  the  stethod  undereatlsiatea  drag  by  a  algnSflcant 
margin  as  sliown  later.  The  cause  can  be  traced.  In  part,  to  the  neglect  of  'higher  order* 
effects  In  thr.  efreomwise  momentum  equation  and  In  the  laetchlng  between  the  vtseoua  and 
Invlscld  flows.  A  revised  version  of  the  program,  known  as  BVOK,  haa  therefore  been  deve* 
loped”  inrl'idifig  tijefle  efferra  together  with  eorrectlona  to  the  lag-entralniaent  isethod 
similar  to  ihosc  in  slVf  described  previously.  (A  slightly  different  shape-parameter 
reianonshlp  from  ihat  of  5ivr  la  used  which  la  considered  to  be  suitable  for  flows  with 

tral 1 Ing-edge  sepup^^lon). 

Drag  Is  calculated  In  BVQK  by  both  the  'local*  and  ’far-fleld*  methods.  However, 
for  r-as»'ns  given  \n  s»ctlon  7,  the  'far-fleld*  method  la  generally  preferred,  and  predic- 
tl'^ns  of  drag  i^y  nvoK  and  VOK  shown  later  have  been  obtained  In  this  way,  using  MACHCOMT 
as  the  subrcMMn<»  fr.r  wave  drag. 

Ftampiec  of  predict  lotis  by  VOK  and  BVhR  of  overall  forces  and  pitching  sioment  are 
sli'’wri  111  Pfg  »  for  n  e*rlws  of  i»i%  thick  aerofolla  with  relatlvely-large  rear  loading. 

Thi"  rig.MA  \n  taVAi,  fpnm  Pef  wh*rw  d»ta1\s  are  glv*n  of  the  aerofoils  and  the  wind 


(•attfunhen  }  (1  ^ip) 


ft)  7  lift  and  (1^*9  versus  inodtort.  HLR  iwo  tlem^ni 
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tunnel  leeaftureiients  u«ed  In  the  taaeaaoient  of  CPD  — .yai 

■ethoda.  Here  It  aufttcea  to  note  that,  at  the  /•r*#**'^* 

lower  of  the  two  chord  Reynolda  nuabera,  /  /*•**  j****^..^-^**-*-  |j 

R  •  6  »  10®,  flow  aeparatlon  la  calculated  by  BVOK  '  f  9  ,f  |  I  S7» 

to  occur  on  the  upper  aurface  of  three  of  the  f  J  \\  [  {  ■  > 

aerofone,  RAE  52?5,  5230  and  523**,  the  ehordwlae  /  /  \V  I  • 

poaltlona  of  the  aeparatlon  point  belne  at  99*,  /  /  •  !• 

95*  and  98*,  reapectlvaly ,  for  C|^  •  0»a.  Hence  »  '  K  ^ — *  (n»V  *  tw  ii^ 

theae  flowa  preaent  a  challenge  to  CFO  aMthoda  for  ^  m  r 

predicting  dreg.  Mr  /  ,/  m  _1'Lr'  //  If! 

Pig  8  reveala  that  the  predictions  of  drag  7  {rrlf  ri 

t>y  BVQK  are  In  good  agreement  with  meeaurement  for  f  //  ({  '  r  li  r< 

flowa  wltti  weak  aiiocka  at  both  Reynolda  number#.  i,  y  xk  \  V  A/  u> 

Therefore,  by  Impllcatlor.  BVQK  predict#  accura-  i  ,y  If  h!  *• 

tely  the  differences  In  drag  between  eeetlona  at  a  *  *  *  *  ^  •  w V»s  » w  • « 

S] ven  Reynolds  numbers  and  between  Reynolds  num-  // 

era  for  a  given  section.  The  Improvement  in  Mr  f  tj 

agreement  with  meaaurement  compared  with  the  pre*  ii  S  •/  II  I 

dlctlona  by  VOK  la  eapaclally  evident  at  •;  jt  t'y  £  P  jiiM* 

R  •  6  •  10^,  where,  aa  noted  before,  eeparatlon  la  as  *  ‘f  (  (  U  f!  I* 

calculated  to  occur  on  the  upper  aurface  of  three  9  ■  _ _  r  i 

of  the  aerofolla.  However,  the  drag  estimatea  by  y  //^  .  .h  f*  !“ 

BVQK  are  less  satlafactory  where  there  la  algnlfl-  ^  ' V  *  •n«h\»  •» 

cant  wave  drag  ( AC(j|f  >  0.001).  Two  possible 

etplanatloiia  are  given  In  Ref  23.  one  related  to  f  r  /  _ n-^-^  /  { 

the  fact  that  HACIlrONT  aaaumea  that  the  local  flow  ’  J  /  ’  |  | 

la  normal  to  the  shock  wave  and  ttie  other  to  the  f  w  fg  /  f  1  f 

tendency  for  BVOK  to  undereatlmate  the  rear  /  /  lx  '  r  I  I 

loading  for  flowa  with  significant  rear  aeparatlon  aiF  /  I  I  i. 

(notably  RAE  5?3<').  'A  study  of  possible  causes  .  J  ,  /  .  .  Lt-u  ‘ 

for  the  latter  ef/ecl  auggests  tijat  the  correctloi»  ^  ••••»#»•  m  no  ui  « 

to  turbulence  structure  for  flow  curvature  is  of  Utt  •••!«»  naoi)S 

doubtful  validity  for  separated  flows  and  la  pro-  a  •  4  •  #  1  ?•  »  w*  hi  l  •  W*  RAt  S2)#i 

bably  best  Ignored  in  such  cases.  The  result  of  neglecting  this  correction  is  sl>cwii  In 
Pig  8  for  RAE  523n,  the  modified  calculation  being  referred  to  aa  -CURV.  Tlie  improved 
predictions  of  rear  loading  wltti  -CURV  lead  to  estimates  of  pitching  moment  and  drag  at 
t))e  'drag  rise*  condition  In  better  agreement  with  measurement. 

A  version  of  VUK  Is  available  with  allowance  for  wing  aweep.  Known  as  SWVGK,  tills 
mettiod^*  represents  the  Influence  of  cross  flow  on  the  shear  layers  but  does  not  include 
errecta  allowed  for  In  BVQK,  whlcli  are  known  to  become  Important  for  unswepl  aerofoils  aa 
separation  la  approached. 

The  accuracy  of  the  predictions  by  thla  method  and  also  by  VOK  and  BVGK  c>r  d.ag 
differences  between  sectloria  and  between  Reynolds  numbers  have  been  studied  by  comparison 
with  data  from  a  pane)  wing  swept  at  25*.  In  thla  aaaeaamant,  the  effect  of  sweep  on  drag 
In  VOK  and  BVGK  la  allowed  In  a  simple  way  aa  discussed  In  Ref  23  which  also  describes  tiie 
aeiofc!)  ae^'t'^na  and  tlie  wind-tunnel  testa.  Here  It  may  be  noted  that  (a)  section  drag 
was  determined  by  the  wake-rake  technique  and  (b)  the  wing  was  cylindrical,  of  sym¬ 
metrical  section  and  waa  t«sted  at  aero  lift* 


lift  Sfii  4*1  siltkmi  •( 
’  Sat  a  #  1  21  a  V  » 
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Comparisons  are  shown  in  Pig  9  between  pre¬ 
dictions  and  measurement  for  the  difference  In  the 
notional  drag-coefficient  per  surface  •  Cp/2 

between  the  two  sections  RAE  5237  and  5238  over  a 
range  of  Mach  numbers.  Theae  sections  are  related 
liirough  calculated  boundary-layer  characterlstlra 
close  to  the  trailing  edge  to  the  unawept  aerofoil 
aectlons  RAE  5225  and  9230  (aee  Pig  8).  Of  the 
ttiree  methods,  the  best  agreement  with  measurement 
la  obtained  with  BVOK,  suggesting  that  the  effects 
shown  to  be  Important  for  unawept  aerofoils  as 
separation  la  approached  have  a  elmllar  signifi¬ 
cance  for  wings  of  moderate  aweep. 

Ttie  effect  on  the  variation  with  Mach 
number  of  the  drag  coefficient  C[j2  of  changing 

chord  Reynolds  numl>er  from  6.5  ■  10^  to  ■  10^ 
la  ahown  In  Pig  10.  Again,  the  closest  estimatea 
of  thla  change  are  obtained  with  BVOK  and  this 
figure  taken  together  with  Pig  6  shows  that  BVOK 
has  a  potentlalTy-uaeful  role  In  the  eitrapolatlon 
of  wind-tunnel  data  to  'full  scale*,  at  least  for 
wlnga  of  moderate  aweep  and  high  aspect  ratio. 

(11)  Methods  baaed  on  the  Euler  equatloni 


A  code  for  the  numerical  solution  of  the 
Euler  equations  baaed  or.  the  flnl te-voluma  method 


....  vcc 
— $wvo« 

- avw 

»o  tlftuwt»l 
f,  a*  («*«er|9l  »»kitwn 

'  asiammit  ><»•  SVvoi 

tkti  nttk  H«  i»i 

RAI  S2)* 

.to,  . 

.  .  i!i- 
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or  et  ^1’'  h^s  h»»n  «rltt«n  at 

Fll^on’*.  To  f'»rm!t  rl^tall^rt  compariaon  with 
a*pwrtm#'nt ,  all'^wane*  haa  been  laade  for  vlacoua 
err»cta  via  the  fneihod  due  to  Wliliama'^ 
le,  using  an  SI  coupling  and  Including  certain 
’hlgh-r-ordcp*  eTreeta,  Drag  la  computed  ualng 
ti»e  Tar-fleld'  method,  the  wave  drag  being 
Inferred  from  th#  loaa  In  total  preaaure  acroaa 
the  ahock  In  th»  way  augg^ated  by  Sella^. 

Only  onmparlaone  with  aieaeurement 

hav»  b»*en  puhllafi^.t  nyf  thene  Indicate  that  the 
laethc'd  gives  accu'a*®  nred  1  ct,  1  ona  of  drag  Tor  the 
sections  RAF  and  o?io  at  high  Reynolds  number 

(Fig  in. 

n#*'AMtly,  Hall’’  has  developed  a  multi-grid 
arhomn  f'^r  snlwliiR  the  F.ulwr  equations  which  whan 
combined  with  ♦e''hnjquea  similar  to  thoae  ater- 
tloned  ahovw  for  solving  the  ahear-Tlow  equations, 
prom1«es  a  meth«^d  for  ca  1 ''u  1  a 1 1  ng  vlacous  tran- 
aonli"  riow«  nvr>r  orrorolls  that  can  represent 
Sh''rk  wav#*  aor-j  r.it  a )  y  whl)#  boing  no  mof*  costly 
t''  t  'v.\  t  ha  I  nv'>  . 
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- SWVCK 
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t  ■»  tatatiM  aklimaWt 
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•t n I iti- t  rvii ’  •  »i»*'  'Ir-.crlh#'!  a  mathod  for 
a''lv1ng  t'l#  n/'yi  i  l---(v#rag**(J,  Navler-Stokca 
r<Tio<’''n''  f 'I  ii^»i*nnn  flow  around  aarofona 

whl'h  la  ha*#’!  -'ll  (he  w''rW  of  W#atlierlll  et  St 
f'-f  fan  1 1  1  p  1  #  acr-^f'-lla.  In  this  method,  Reynolds 
f'leaawfi  ar#  lining  Hi*  *ddy-v  Isco*  1  ty 

I  ly  t  1  ho  a  1  a  rnfril' 1  ii«*d  wpit  an  alg*  brolc  turbulence- 
ir-  1'').  Ti)u»  »h#  Ta<~»!i--d  la  probably  not  rallabi# 
T'-t  pr^  ’t-tlng  -’i  Ip  r  >r  cn*'**?’  with  ragtona  of 
*  nr  ^  1  "U  unar  Mm  trailing  #dg#  au'^h  as  thoea 
.«  -sfm  1  r  -1  I'rn','  J  -  1  V  . 
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Iff  nr  M  - c-r  wlng'^  ar#  aithar  Invlacld  or 
ain  /-r  Mro  vi-r-n-'**  invlscld  lntara<*L1on  type.  Th* 
vIcnMin  f  thaa*  m®fh«da  ara  HOt  SS 

atj  .^i'cad  ai  Mi''«n  f  r.r  saf  rfAiia  In  the  treatment 
of  •ffarin  wM  '  h  ar#  «lgMiri'*ant  for  fJowS  thSt 
ara  cloaa  t anp^raMou  and  ''onaeqijant  1  y  cannot 
yat  pradl'i  rl|0  -Itag  r,f  nodarn  wlngs  wl»h  the 
firritrfiry  d*m/>.jin t  ra t o.i  in  Figs  8  to  11.  Generally, 
tl,r.  vinnn,,,  „  r  \  .|aa  I'lraCt  COUpllng,  although 
01  no.jpjif.p  t«  nn{l/'yad  In  su  appooilmala  way  In 
nne  nnM-.r,d  (aan  latar).  Daaptta  lacking  the  acct|- 
racy  o(  tt,#  matitoda,  wing  tachnlques, 

ii«»d  wl*»,  .-a  it  Ire  and  wfpnrlanca,  ara  Invsluible 
si'lr  •ir'-igi'.  fr-vl'llng  th#  facility  to  Identify 
Olid  niiiliflirrt  1 1,  f  or..  ,1 1  ni#iia  1  ona  1  aonrcoa  of  e*cess 
rlr  ar . 


•"Ig  1’  1 1 1  u  I  *>  ♦  «*>  tiio  mathods.  or  the  pan»l 
.1^  7^,  »' .  1  I  .  »  J  ♦!,«»  ,1,,*  fatrle*?  (.’?rARV) 
app'’ara  [•'  br»  n,n  rri''**  uao’l  Slid  Is  th*  SUbJeCt  Of 
co/»r  1 1)./ 1  rip  I ,  Allowancs  is  Included  In 

till*  mofit^i  frr  th'’  afTort  of  wing  boundary  layers". 
Tho  ffivf*r*ii  tren^oiiir,  ema  1  I -part  u  rba  1 1  on  aiethod 
r-f  aiboiir’  at  « 1  '  '  witli  vlscoya  affects  Incor- 

porafad  <’y  Flrmlii",  |s  now  Mrgajy  suparsaded  by 

tha  mora  ac'*ufaie  f  u  I  I  -  pot  ant  1  a  1  and  Euler 
waihod".  Tiia  f  11 1 1 -potaiif  la  1  mathod  of  Forgay  and 
rarr**  (Ff)  has  haa,i  usnd  for  savaril  yatrs  and  \a 
gonarally  t-agardod  ea  a  good  axsmpla  of  a  method 
pf  this  tvpa.  A  version  of  the  aiethod,  due  to 
Arthur,”  la  avalisMa  with  allowanca  for  viscous 
errerts  (VFp).  Finally,  BAa  Fllton  hava 
programmad  a  t  h  r'aa- ill  mans  1  ona  1  varalon  of  the 
F.oJar  aiaihod  raferrad  to  In  section  tn  thla 

method’^,  tiia  sh-ar  layars  ar*  calculated  on  the 
asaump»jr*n  of  plsfiar  flow  at  each  streamwlse  aec- 
1 1  on  witi.  tiio  saiiitinn  coiipled  to  the  Invlacld- 
rir>w  sr.*v-i'  bv  an  f  schamw. 


■  A|  »» 

ft  r  ?S  • 


fi|ii  vsrteima  ftt  Sr»|  itftfaw^t  wiik  an 
fftfftaitat.  n  •  1 7?s 


Fow  raS' 
a  t|-n  r  1  mr*  n*  s  ra 
r  nnr*#  r  na  l  w  M  li 


iM*  nf  compartaons  of  drag  predictions  by  theae  laethrda  with  wind-tunnel 
nvsiiai'la  I'nr  puhllcatlnn,  and  consequently  the  remainder  of  the  section 
ffioMw.d*  of  analyalng  the  drag  of  wings  from  Intormstlon  provided  by  the 
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codes  based  on  ttie  cleaslcal  Tar-^fleld'  approach  described  In  section  2.  Results  or 
lysis  are  presented  to  Jllustrste  the  power  of  this  approach  in  identifying  sources  of 
excess  drag. 


iJ  >1—1 


An  analysis  of  drag  la  shown  for  a 
wing/body  suitable  for  a  transport  aircraft 
comprising  a  wing  of  aspect  ratio  8.  with  a 
leading-edge  sweep  of  2e*  and  a  trailing  edge 
sweep  outboard  of  the  tralllng-edga  crank  of  14* 
(Pig  )3)>  In  addition,  a  study  of  wave  drag  la 
presented  for  a  wing  representative  of  that  of  a 
subsonic  combat-aircraft  having  leading  and 
t'*9lHng  edge  sweeps  39*  and  15*  and  an  aspect 
ratio  3.3- 


(1) 


Transport  aircraft  configuration 


Comprehensive  CPD  calculations  are  not 
available  for  this  confieuratlon  and  so  the  analy¬ 
sts  is  performed  using  wing  surface  pressures 
measured  on  a  complete  model*.  Limited  ealcula- 
tlons  of  wine  pressures  for  this  configuration  by 
both  the  RAe'*  and  VPP^^  codes  have  been  tound  to 
be  In  reasonable  agreement  with  sieaeuresient  (made 
In  the  latter  case  on  a  related  half  model). 

The  Vorm  of  analysis  is  Illustrated  In 
Pig  14.  The  body-drag  coefficient  C^g  la  deter¬ 
mined  from  tests  on  the  body  alone,  thereby 
avoiding  the  difflrtHty  of  determining  sting 
interference.  Mole  that.  In  choosing  tlte  ordinate 
for  ttila  figure,  use  is  made  of  tiie  fact  that  tiie 
vortex  drag  is  close  to  tlie  minimum  value  for  a 
planar  wing  »>y  sut-tractlng  from  the  drag  coef¬ 
ficient  ,  wtiere  ^  Is  wing  aspect  ratio 

and  auffl*  PAl.  refera  lu  force- 
bala<..j  mea  su  I  e'i<ent  .  Tiie  small 
exctsa  vortex-drag  c-ieff  1  o t ent  . 

ACoTV  •  ‘^DTV 

la  determined  from  tiie  measured 
apan  loadings  using  the  classical 
Trefft-plane  meltiod  referred  to  In 
section  2,  Two  alternative 
vortex-trace  mo<iel8  have  been  con- 
eldered.  one  allowing  f-..-  the  body 
in  a  etmple  way  and  the  ottier 
representing  the  trice  aa  a  planar 
slit  of  the  same  span  ae  the  wing. 

The  latter  model  was  chosen  for 
the  analysis  on  the  basis  that  It 
yields  values  of  overall  lift  in 
closer  agreeincnt  with  the  balance- 
messur^d  values  than  those  of  the 
other  itodel.  However,  the  excess  vortex  drags 
given  by  tlie  two  models  do  not  differ  by  much 
<  0.0002)  suggesting  that,  where  overall 
lift  ia  known  accurately  from  aome  other  aource 
(In  this  case  the  force  balance),  the  drag 
stialyala  is  not  aensltlve  to  the  ahape  of  the 
vortex  trace. 

Calculated  values  of  sCq’I'v  shown  in 

Pig  15  plotted  against  lift  coefficient  for 
varloua  Mach  numbera.  Except  where  there  ta  s 
rapid  increaae  In  vortex  drag  with  lift,  the 
excess  vortex  drag  varies  alowly  with  botli  lift 
and  Mach  number,  the  sudden  Increase  being  attri¬ 
buted  to  the  ioaa  In  lift  on  part  of  tha  outar 
wing  following  flow  braakdown. 

Except  in  special  eaaas,  tha  Intagrand  of 
the  vortex-drag  Integral  or  'local*  vortex  drag 
cannot  be  related  to  sectional  drag;  however  there 
la  a  direct  relationship  between  'local'  vortex 
drag  and  span  loading,  and.  In  the  present  cate, 
the  cause  of  the  non-zero  excess  vortex-drag  la 
that  the  outer  wing  la  relatively  lightly  loaded 
compared  with  the  Ideal  elliptic  loading. 
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kn  In  w^ll  Wimwi),  th#  vor¬ 
tex  or  wlnR9  with  non-plan%r 

vorto*  tra**^^  (*6  wtng/w|npl^t 
conTJRurat  lo»>9 )  c^n  h^low  tho 
mlnlwu"  for  planar  wlnae  oT  the 
aama  arf^n,  anrt  a  technique  for 
calculat’lng  the  mtnlmum  vortex 
drag  of  non-planar  confl j:\iratlon 
hao  ha»»n  proRrai«m*'d  by  laaaea'*. 


kn  In  ^♦‘rtlnn  ?,  ,  "77^ — 

vlaron*»  drap.  la  lt»(*»rred  from  * 

bouM'la ry- 1  ay**r  <^uan»<t1aa  at  the  *  / 

tratMnR  rdija  ualna  ati  a*tand^d  ^ 

version  of  tha  "-pM  re-Yontig 

Tormnla*.  Th»*  tvnhiil»nt  h^ijndary-  L— _ _  _  _  _  _ . 

layers  ar*  cal'''ilated  using  the  •  .  *'  V  V  *  \ 

eiaaanr^d  rr^aanre  rtlatrlhuMona  ®  T  *  •  •  •  *• 

«n-1  .n  'inrinltr  fr-red  .Ing'  '■*« 

veralnn'^  of  the  1 ag>ent ra 1 nment  method.  Comparlaona  mlth  the  potentially,  more-accurate, 
thre^-dtm-nslnnal  oP  Smith'*"  auggesta  that  the  ‘Infinite  tapered-«ing*  method  almulatee 
n»^r..^Ma •■*'1  y  thrr'o- <n in«na  1  nna  1  •ffeota  In  the  present  caae  except  close  to  the  tip  and  the 

r'^nf  . 


TvplraJ  apanwlae  (^la^^1hu- 

tl^r»r.  r\l'  vl  a»“rnr-drap.  roef- 

rirlon'  lllua^^a^«'1  1  ti 

P1p  )f',  t-'i-  II  '  n,7fl.  Th'- 

r  •  I  a  M  v'' 1  y  •>  I  n  rge  jurf^ari*  In  lAr*al 
V  1  •  r  on  *»  -  dra  p  m*»rt  <''*ent  on  t  !>• 
m^er-  wlnr  lo  in^  oo»fr<clent 
from  0.11,“  tn  n ,  r,^, 
ron«l*»t*»nt  wP'i  tl|e  grOaM>  J  II 
«ih-Nrv  aire>\p«li  wlMi  IIT^  and  M»e 
rnnso'jiiont  tliloU^nltip  of  tli® 

♦  •«ijnd^rv  lay^r  !■  yii“'  i^am  of  th^ 
ahi^rh  on  M»!''  pai*  o(  tti«  wlnp. 
Tno  mapin  M  idr>«  of  Min  lora) 
ornil  rl  t'nM  oii«  t nv«»r»ll  vlaron* 
di  ap  a*o  1  ndl  ' "d  1*’  flc  by 
r  f  .d  o  (  n  V  o  .  wlmt-n  c  1«  lo-'a] 
«||^a<nw]«n  nhnrd  and  <■  I 
pnnifin  t  r  1  n  moan  r  ln>  i  d  . 


In  Ih*  sO'ion'*''  '^r  rio«fi«i«i 
1  nfoirra  wa/o  diap  haa  been 

oalriiiatod  oy  loou'a  method'*’.  It  will  be  recalled  from  section  ?,  that,  In  this  method, 
to*  varjaM'^n  of  ni.orv  atrengMi  with  distance  normal  to  the  wing  surface  Is  determined  by 
wing  atreamwiaa  oijrvatura  end  static  pressure  at  a  point  Just  upstream  of  the  ahock.  This 
!•  to  Ignoring  the  •ff^ct  on  flow  curvature  of  the  boundary  layer  and  assuming 

that  th*  airenr'h  the  rhoob  in  the  field  is  unaffected  by  the  variation  of  surface  cur- 
vainre  -iioijg  tho  oh-^rd  tipatf-eam  of  the  Shock.  These  aspects  are  considered  again  In  the 
fiooond  o?ain| do  fn  wh|oh  there  ift  a  rapid  variation  of  streasMlae  curvature  ahead  of  the 
fiho/i'  on  p-^rt  of  I  ho  win*.  Howover,  In  the  present  case,  the  curvature  of  each  wlaig  sec- 
t|..ii  ie  olnoo  ♦«  a  minimum  In  the  region  of  the  ahock. 


**.pshwj-o  v«r1»tionr  of  the  local  wave-drag 
ooorriofent  o*|o.,]af«d  py  Lock's  method  are 

sii-wii  1,.  Fi*  17  togothor  wlth^the  local  contrlbu- 
t|r»rt  I"  wave  flfor  rpj^f»>)o(»))/c  for  M  ■  0,7B,  The 
cent  rlli'jt  Ion  to  wnv*  dras  of  the  part  of  the  win® 
fnb^wrd  0^  the  r r» 1  1 1 n*>»dy»  crank  la  seen  to  be 
ro)»Mv«1y  wrmM.  with  mowt  of  the  wave  dra«  orl- 
pln^ffnp  ffom  n  fr>g_\nn  Jnwt  o'ltbeard  of  the  crank. 

noth  loos  I  viscous  and  wave  drags  have  been 
Ifiirgrated  serosa  Ihe  wing  span  and  have  then  been 
rnmi.fnod  wfth  vortwi  drag  and  body  drsg  as  shown 
In  Fig  I 't  to  give  overall  drag.  Comparisons  bet- 
wewn  'ra  lei|  fated*  and  mwaanred  overall  drags  are 
shown  In  FIs  and  Indicate  that,  for  aubcrltlcal 
flows  oi-  t»t  Ihe  rogfon  of  minimum  drag,  the 
* ca  1  rii  1  at e.i •  drag  roerrirlen*  1*  lower  tt»an  the 
measured  vain*  I'y  an  amount  which  varies  between 
0.nnn»i  nt  M  -  h.^  and  n.oOOfi  at  H  •  0.8,  aithou|h 
fn  lose  *<•’■•*1  wsieemefit  with  measurement  than  BVO* 
Is  foujid  to  bs  for  a  series  of  aerofoil"  (Fig  8), 
^f.tlmsfos  are  enoonraglngly  rJose  to 
messiirement  and  siiow  that  the  'far-fleld'  aiethod 
hai  •  user,,!  r^'le  M  play  In  tlie  analyala  of  drag 
of  wlng^'hody  oonflgurations  suitable  for  transport 
alrersit.  A  study  the  sources  of  the 


fif  II  Sftswnt  anttbvhss  st  ammtktti  Istsl 
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dj«crep«incl«»"  that  the  error#  can  be 

largely  e*pla1ne<1  hy  Dow  Teaturee  not  represented 
minpi^t^ly  in  th*  Mnalyaja  including: 

in)  wlng/bn»l.Vt  tmnndary-layer  Interference; 

<h)  fi»^w  curvature  an<l  Reynolds  normal  atreaaea 

tti  the  turhulent  shear  layers;  and 

(r)  tranaltion-trlp  drag. 

Pig  IR  reveals  that  the  dirfereneea  between 
'calrulated*  s?id  in««aaured  drags  decrease  as  wave 
drag  Inrresses  frir  Mach  tmaihers  In  the  range  0,7 
to  0.6l.  Tlie  moat  likely  explanation  for  this  Is 
tiiat  f.ock's  methofi  overestimates  wave  drag»  since 
It  la  unlikely  tlial  the  estlmstes  of  the  other  two 
drag  components  become  store  sccurate  as  shock 
strength  Increases.  On  the  evidence  of  studies  of 
jnviscld.  two  dimensional  flows  It  is  stated  In 
Ref  Ro  thal  estimates  by  Lock*a  method  are  pro¬ 
bably  within  -10  to  30J  of  the  correct  value 
except  at  lew  vsluea  of  C()^#  (,  *  0.0015)  when  It 
ceuld  be  up  to  0,no05  too  high.  No  direct  evi¬ 
dence  Is  available  on  the  effects  of  the  boundary 
layer  or  three-d 1 mens  1 ons 1 1 1 les  In  the  flow. 
Mowevet ,  some  compartaons  have  been  made  between 
pr»*'l  1  et  1  ons  by  l.ork*s  method  and  thoae  of 
ajlwrlght's  field  method*.  In  each  case  based  on 
calculations  by  the  FT  method  of  Forsey  and  Carr* 
for  the  present  eonf 1 gurst 1  on .  These  comparisons 
reveal  that  threr.-d  1  mens  1  ons  I  effects  are  slgnlft- 
csni  only  In  tlie  ne.sr  vicinity  of  tt>e  tip  (le 
wivhin  shout  r»ne  nr  two  chords)  and  thus,  overall, 
thoii-  inriitonra  ou  wTvo  drag  may  he  Ignored. 

(ID  Combat  a  I  rcraft  wing 


MIHdlMIN* 


•  0.*5  »  Lock's  method 


Ttie  second  ronflgiiratton  la  an  example  of  s  wing  design  for  which  Lock's  s»ethod  - 
nt  loftsi  In  its  present  form  -  is  less  reliable.  The  wing  has  been  tested  as  a  half  mode) 
wM»,  the  Sim  of  providing  fl u  1  d-dynsmic  data  for  the  validation  of  CPD  methods, 
cotrtpai  1  sons  ef  predjotions  by  VFP  and  measurements  of  wing  pressure  distributions  are 
dis.iisscd  tn  Ref  q.  As  part  of  this  study.  N.  C.  P.  Flrmln  (RAF)  has  performed  some 
ca  1 '■'I I Ions  of  wsve  drag  using  lioth  Lock's  and  Allwright's  techniques.  Results  for  local 
wave. drag  aie  shown  tn  Ftg  IQ.  f>»jtbosrd  of  the  shock  bifurcation  at 
is  seep  (o  give  imjrh  larger  values  of  local  wave- 
dreg  than  those  of  Allwrlght  while,  further 
Inboard,  Lock's  predictions  for  the  resr  shock  are 
slightly  lower  on  average  than  Allwright's  values. 

An  explanation  for  the  former  discrepancy  la  given 
In  Fig  pn  which  ahowa  the  variation  with  distance 
from  and  normal  to  the  wing  surface  of  shock- 
upstream  Maoh  numher. 


A I  s  »  O.hn'i  ,  ie  outboard  of  the  bifur¬ 
cation,  lock’s  method  predict#  that  the  shock 
penetrates  m*M‘h  rurther  Into  the  field  than  Is 
Indicated  |iy  the  more.a Ceil rs t •  field  s»ethod  of 
Allwrlehl.  The  reason  for  this  te  that  the  cur¬ 
vature  of  the  wing  urper-svirfsce  Increases 
ms»kedly  with  dDtane*  upstresm  of  the  shock  on 
this  pof I  of  the  wing.  Thus  the  flow  curvature  at 
the  shocir  in  the  field  is  affected  (vls  the 
outgoing  Alacft  eharacterlstics  from  the  wing  aor- 
fsre)  and  c'^nsequent  1  y  the  rat*  at  wtileh  Rji 
changes  wlih  dintaiice  nnrmal  to  the  wing  la 
modified. 


Fig  ?n  also  shows  that,  close  to  the  wing 
eurfsce,  wliere  the  flow  la  strongly  Influenced  by 
conditions  at  the  foot  of  the  atiock,  there  la  a 
marked  difference  in  the  two  predictions  of  the 
variation  of  "H  with  distance  from  the  wing. 

This  discrepancy  arises  from  tt»e  neglect  of  the 
effect  of  tlie  boundary  layer  on  (a)  the  local  flow 
curvature  and  (h)  the  inclination  of  the  shock 
relative  to  the  wing  surface. 
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Oewptt*  thpfl**  ft<»ric1»Mc!^8,  Lock*8  Method 
\n  Meerul  In  provHlng;  «  r^pM  Indication  of  80ur*> 
ce"  of  e*c**«8  drag  both  In  th*  atagea  of  the 

wing  denlgn  and  on  aa  a  dtagnoatlc  tool 

following  wlnd-tunnol  twata. 

3 . 3  Bodlea 

ferhaps  tl»o  flrat  UK  attempt  to  uae.CPD  for 
tho  prediction  of  body  drag  waa  by  Hyrlng  who 
employed  a  viecona/lnwlacld  interaction  technique 
to  cal'*«ilate  the  anborit.\r.al  flow  over  axtayM' 
Metric  bodlee  7ero  Incidence.  He  repreaented 
the  ihvlacld  rinw  over  the  diaplaceMent  surface  of 
the  body  and  th«  .che.ir  layer  by  a  aource^rlhg 
MetJi^d  and  calculated  the  vlacous  ahear-layera  by 
Integral  meth^de,  coupling  the  two  aolutlona  by  a 
Plrect.  procedure. 

UaJng  hla  method,  Myrlng  waa  able  to 
deaian  a  Mow-drag’  body,  aa  Illustrated  In 
Pig  ?1,  where  tt  la  dl at  1 ngulahed  from  a  conven¬ 
tional  body  of  th«  same  thlcVnesa  ratio  in  having 
nn  pfVHncuMce/]  etjcijon  peake.  Alao  shown  In  this 
figure  la  variation  wltlt  tIHckneaa  ratio  of 

drag-cnei  r  1  cirtif  hna«d  on  aurfice  area,  C[)|^  ,  for 
both  typcf*  of  body,  cJ»ai-ly  tlluatratlng  the 
auperlorlty  or  the  Mow  drag’  dealgn,  albeit  at 
the  r»r  a  low«r  body- vol  umc .  On  the  other 

hand,  »hc  'low-drop'  body  hap  somewhat  higher 
aucitona  or  local  veincntea  than  those  of  the 
cc-uyeMMonnl  ehcp/>  \u  the  region  where  the  wlnga 
of  a»)  aircraft  might  be  mounted,  showing  the 
danger  r>r  op »  ! m j f nP!  aircraft  componenta  In 

1  -jel  a  »  1  nn  , 


A  nnmbc,  mcH,nda  have  been  developed  In 
UV.  fop  ra  I  cij  1  •»>,  t  ng  irat^aonlc  flows  over  bodies, 
tnc}tidJ{jg  »  he  }  u  I  t -pot  «*nr  1  al  method  of  Baker  and 
ngie’*'  lor  »x  1  aynmot »  1  p  bodlea  and  two  methods  of 
anivliig  the  Tulor  erpiafiojia  for  the  flow  over 
f  c  r^bnd  1  cr<  . 

hakcr’a  method  \\nn  been  used’’*’  to  calculate 
tho  vai’1at\on  of  drag  with  Mach  nnmher  of 
aphertca  I  t,v-bi«jnted  fopebodlca  Bt  zero  Incidence 
for  Mach  uumbera  tip  to  the  limit  of  validity  of.^ 
the  method,  le  app I'o *  1  ma t e  1  y  Unity.  An  example 
of  the  reasonable  agreement  between  predictions  by 
this  method  of  presbure  d 1 e I rl but Iona  And  drag  la 
provided  by  pi^,  (irag  la  Inferred  from  the 

calcijij^tirtM  by  the  Moca)'  method  and  a  small  tare 
correottou  to  a)|ow  for  discretisation  errors  In 
tl»e  methofi  -^nd  ab  I II- f  r  i  C  t  i  Oft  dreg  Is  applied  to 
tbe  thenry  to  align  prediction  and  measurement  at 
M  -  ».7. 


•  I  f7  •«  **f/(*^ 


It  CakaUlMa  *1  Hxiti  itvwtti  nmaM  It  shteX 
tl  let  tfthtnt  ta  wM| 


fl|  II  tithitvhent  tnt  tTt|  itrtm 

tMtlinfff  r«h0  fm  tclifwwrhli  VvStct 
<*lt»l|l»t  tf  Hfrlnit  wtllitt 


forr^ap-^ndlng  calculations  of  drag  by  the 
rir«t  of  the  F.ulor  methods'*^  are  also  shown  In 
Fig  ?7.  Paeed  on  I  lie  BAe  algorithm  fOP  aolvlng 
the  RuJep  equntiona,  this  technique  la  applicable 
to  aileymmefrlc  f opohodl es .  Again  the  predicted 
varldfloii  of  drag  (by  the  Mocal’  method)  with 
Mn^'h  nnmhcr  In  the  nubannlc  range  la  In  fair 
agreement  with  measurement.  This  method  haa  been 
gcnepaJ1ao»i  by  RAe^*  to  Include  forebodlea  of 

fenepsi  ehape  at  Incidence,  and  a  further  genera- 
laatjon  hao  been  fiepformed  hy  Aircraft  Research 
Assoclatlnn,  (ADA)  Retlrord***  who  have  applied  their 
Mulflbl'^rk  fechulque  to  enable  sideslip  to  be  con¬ 
sidered.  freanure  distributions  on  the  upper  and 
lower  alden  of  the  ^dy  calculated  by  the  latter 
Method  ape  compared  with  s»eaaureMent  for  the  fore- 
body  of  the  RAe  Hawk  at  incidence  and  aldesllr 
angle  g  In  Fig  ?3*  No  eompariaona  of  drag  are 
available  hut  ttic  agreement  between  calculated  and 
Measured  pressure  di n t r 1  but i ons  la  reifonablr  good, 
Suggwatlng  I  hat  tlic  method  may  be  used  to  calcu¬ 
late  tfie  vartaMoii  of  rfrag  with  Mach  number  for 

auc|»  nhapcc  . 


Sm*'.  tm  rn***’** 
ikt  ie*» 

»  iim«ww>» _ 
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*eyi  XffMwt  eii'e«>w«i  ww  e«f  r«»  •  rMaw 
tfvHTM  I  ipat'iiii  «•«!  Ww»iw«  «iiw 
r/a  >  «| 


15 


Techn)q<i*A  such  an  the  last  one  have  ]ret  to 
be  c<»»Mne<J  Interactively  with  boundary-layer 
calculation  methoda  to  predict  the  drag  of  general 
bodtea.  Of  partle«ilar  Interest  in  this  connection 
are  r»eei*gea  wtth  upswept  sTterbodlea, 

3 •  Cowls  _ajid  .n<Mi_2 l^a 

The  accurate  ralculatlon  or  turbofan  cowl 
drsK  Is  an  Importaut  rone iderat.  1  on  In  the  design 
and  the  perrermance  predlettnn  of  laodern  transport- 
alrcrsrt.  To  he  fully  representative,  the  calcu¬ 
lation  method  should  simulate  the  Interaetlon 
between  the  engine,  the  pylon  and  the  wing.  This 
cannot  done,  st  present,  although  progress  la 
being  made  in  the  s»odell1ng  of  complex  con¬ 
figurations  (section  3.0)  but,  aa  a  preliminary  to 
obtaining  solution  to  the  complete  problem,  two 
Biethodn  have  hoen  programmed  for  Isolated  cowls. 
These  methods  have  a  similar  function  to  that  of 
aerofoil  methods  in  providing  a  simple  basie  for 
checking  flow  algorithms.  The  first  method,  due 
to  rears'*’,  us*s  a  Plrect  coupling  of  a  full- 
potential  solution  of  the  Invlscld  flow  with  the 
Isg-ent ralnment  isethod  for  the  turbulent  boundary 
layers.  The  aecond  procedure  replaces  Peace's 
potential-flow  scheme  by  the  PAe  method  for 
solvinc  the  equst Ions ^ 

Ooldsfni  th**  •  liso  made  compsrlaona  between 
rt  1  rxnf  hy  these  methods  snd  mess»irementa  of 
cowl  pfess<if'«  ftrsg  for  s  numh«r  of  NACA-1  cowls 
aligned  with  the  ft-nr  -tresm.  Comparisons  for  the 
cowl  ge-Mnrti'v  sket  clio.j  \i\  ?*l  sre  Shown  In 

Pig.  where  rowl  f'ressufw-drsg  coerflclent  la 
plot*e,i  Agalnni  the  reMt  I  ve-flow  ratio  aa 

defiii*i.i  111  Fig  ?N.  rwflce's  siettiod  la  limited  to 
Hseh  iiMmhers  he|ow  shout  of  unity,  and  In  this 
Msch-niimher  range  it.  gives  good  sgreement  with 
m«s«'irement  Tor  r*  1  s t  f  ve- f  I ow  ratios  shove  those 
for  whief,  cowl-lip  sepnratlon  occurs.  For  low 
re)s  M  ve-r  1 ''w  rotlos,  the  agreement  la  less 
sa 1 1  srs'^t nry ,  as  might  he  expected  for  a  method 
using  a  first-order  treatment  of  the  shear  layers. 

The  pMler  method  has  only  been  used  for 
ca  1  eiii  at  1  ons  at  supersonic  speeds  and  so  a 
dlscosslen  (,r  itiesc  comparlsons  Is  deferred  until 
section  S  where  msthnds  for  supersonic  flows  are 
discussed. 


ri|ll  Vxatr  Ss4  hvir  Ml  r*««vre  g|(ra«htn». 
lAi  Hawk  ffttfbtSy.  H  •  l.|,  ■  ■  1.77*. 

D  •  tii* 


*1  '  f*wl  er  tft 

•  «»,»  1  .  kU  A.i,  » 

,, 

*  7“  [  Is'*' 

J 

u 


ti|7i  C»vl  |»»wvlrr  intf  rtthnttlent 


A  number  of  stethods  have 
been  produced  In  UK  to  calculate 
the  di-ag  of  afterbodies  with  pro¬ 
pulsive  Jets,  Hodges'*’  has  con¬ 
sidered  111*  esse  of  an 
sxisymmelilr  aftertuMly  with  a 
single  fni  Slid  simulates  thr 
external  flow  hy  s  panel  method, 
the  Jet  hy  the  method  of  ctisrao- 
terlstlcs  and  I  he  boundary  layer 
with  the  Isg-entrsinment  stethod. 
Thus  the  mol  hod  Is  restricted  to 
uni f rrmi y -subson I e  external  flows 
and  Jet  flows  which  are  entirely 
supersonle.  The  solutions  to  the 
various  parts  of  H»e  flowfleld  are 
patched  and  emplrlea!  relationships 
are  used  to  define  the  separation 
and  resttsrhment  points  snd  also 
the  entrainment  In  the  mixing 
region.  Comparlcnns  of  prediction 
by  two  stethods'”,  Including 
Hodges'  aiethoil,  and  sieasuresients 
of  afterbody  pressure  drag  for  a 
series  of  noxsJes,  at  various  Jet- 
prasanre  ratios  end  for  M  »  0.6 
and  O.A,  reveal  that  Hodges' 
method  la  In  ressonsMe  agreement 
wtth  messiirement  for  sul'Cr*  t  Ics  1 
external  flowa. 
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ha^  d^v^loped  a  mathod  baaed  on  aolutlona  of  the  Euler  equatione  In  both 
th»  external  flow  and  the  Jet  which  la  not  reatrtcted  to  suberltleal  flowa  outside 
the  Jet.  As  In  MAdgea*  procedure,  the  boundary  layer  la  calculated  by  the  lag-entraln«ent 
method* >  hut  rerlanes  the  Olrect  coupling  and  empirical  aeparatlon  prediction  of  Hodgea' 
technlo«ie  hy  an  ??l  coupling.  On  the  otJier  hand,  the  entrainment  In  the  Jet  mixing  region 
la  determined  hy  a  afmple  empirical  correlation. 


P|g  shows  plots  or  afterbody  pressure- 
drag  coerriclent  against  free-atream  Hach  number 
for  an  afterbody  nottle  configuration  tested  hy 
neubueh  and  nunrkei'>.  Tlje  predlctlone  by  the 
meth-^d  nf  Peace  ere  neen  Vo  be  In  good  agreem*nt- 
wUh  meaniiremeuf  eitrnpt  rjoae  to  H  •>  1 . 

3 . S  Complex  conf 1 gu rat Iona 

Th"  reqt»l  remn,tt  to  he  able  to  calculate 
transonic  flows  around  complex  ro»«f  Igura  1 1  ons  , 
such  as  those  a»»own  in  Pig  ?7,  has  led  to  the 
development  at  APA,  P-df ord' ^ and  at  EAe  of 
•ultIbJork  grid  generation  schemes.  Combined 
with  the  pxe  technique  for  solvlfig  the  Euler 
equations,  these  methods  lisv®  been  used  for  the 
calruJstinn  of  the  riow  over  a  wide  variety  of 
conf  Iguj'st  1 'uip  ,  an  ej-^^mple  being  given  In 
section  'I.  ffnwever,  assessment  of  drag  predic¬ 
tions  by  the  method  Is  still  at  an  early  stage, 
and,  as  ji'^Ved  in  se'*tlon  ?,  the  production  oi 
spurl'Mis  etitropy  hy  M,s  current  generation  of 
Euler  solvers  makes  Ihe  accurate  determination  of 
drag  dlfflculi;  nevertheless  it  is  envisaged  that 
p'Nsslble  applirst  of  the  method  In  the  future 

1  nr  I ude ; 


(1)  d"  t  m  1  iia  1 1  fMi  of  the  Installed  drag  of 

Pvi<-M’>cv(i  r>r  wr.np«n  srraugements ; 

Mil  r>f  tnnio^l  drag  cf  closely- 

»•  «Mr  1 1  '•'•'uf  ^  p.M  fs  t  i  ;  and 

(Ml)  ra1c,||sM''M  drag  wlng/wlnglet 
oniM  M*  ♦  1  . 

A  MF.ninhr!  pun  -.iirFn.cotHr  AinrnAPT 


The  s'rframe  oAmponents  of  supersonic 
aircraft  ar*  generally  Integrated  closely  and 
hepoe  the  rody nam 1 1"  InterTwrence  between  them 

can  be  coftslderahlo.  ronseq  jently  this  section  is 
different  from  the  f'receding  section  In  that  no 
d I s 1 1  O'" t  1  ou  Is  ma<le  br>tweeu  compofjents  and  the 
mei  hods  a»-e  con«idor««l  under  separate  headings  in 

rhroir^l  r«p;i --1^  1  ^rder  np  development. 


Generalised  near  field  wave  drag 


The  discovery  that  met.hods  based  on  'area  1*^7^  Swlact  9fal<  nt 

thsnsfer’  rules  do  »>ot  glv*  reliable  predictions  tonltaur^hons 

of  x»ro-))rt  wave  drag  led  BAe  (Warton)  to  produce 
a  codo  has»d  on  a  simplified  panel 
method  for  Ilpenrlecd  siipersonlo- 

flow  known  as  the  r,oi,*ra  1 1  aed  Near  Pssstsfs  tNp<nf  Is  Ws^cs 

FlPlJ  Wnv.  Urn.  ('IIIFWIM  prfFr.m".  li».rwwlc  Iw  tW  Bn,  (C«J 

•^nffirlenl  oonridourr  J)*r.  hecft  -  ■  Cmpv 


eatat’llahed  in  the  accuracy  of  the 
method  for  a  range  of  military 
combat-aircraft  configurations  for 
It.  to  he  used  In  a  routine  way  ot> 
pro.lect  design.  An  appMcatlon  la 
Illustrated  In  Fig  the  design 

sxerclae  Involved  changing  the 
fuselage  geometry  and  estimating 
drag  ualng  the  procedure.  The 
particular  dealg»»  alteration  ahown 
In  Pig  P8  Incresaed  fuselagf 
voltimo  while  redurtng  aero-ifft 
drag  I’y  IJ*.  A  comhtnat.lon  of 
changers,  auoh  as  at  ra  1  ght «til ng  the 
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8pln»,  wAlntins  ih*-  nlfleSt  8nd  Increasing  centre- fuselage  volume,  reduced  xero 

lift  dr8g  hy  sf  end  lT>rr»«8*vd  Internal  fuselage  volume  for  fuel  system  etc  hy  HoO  litres. 

/tl though  largely  superseded  by  recent  developments  In  methods  for  solving  the  Euler 
equations,  terhnlquen  such  as  ONPWD  currently  retain  an  Important  function  In  the  design 
of  supersonic  combat-aircraft  because  they  are 

(a)  economical  In  terms  of  computer  time  and 

(b)  simp]*  to  MS*  and  to  understand. 

*1.?  Eule_r/J>anel  program  for  wlna/body  eonf leuratlona 

WhlJ*  giving  r*asnnah1*  predictions  for  the  flow  over  bodies,  panel  methods  are  not 
aatlafactory  for  lifting  sucfnces.  In  general.  Thua  a  hybrid  procedure  has  been 
programmed  by  PA*  (Warton)  using  an  Euler  code  for  the  wing  and  a  panel  method  for  the 
body55.  The  method  haa  been  used  to  predict  incremental  drags  and  pitching  momenta  (from 
eurfac*  pressure  Integrations)  for  a  combat  aircraft  configuration  due  to  wing  camber  and 
twist.  Pig  PO  shown  that  the  procedure  gives  accurate  predictions  of  the  changes  between 
two  rtirr*reMt  winga  over  a  relevant  range  of  lift  eoeff iclenta . 

Trim  drag  In  an  Impnrtant.  rnna  1  d*rat Ion  In  the  design  of 
combat  airrraft  for  both  subsonic  and  supersonic 
fflario*ovr*a.  ^*oMssnu*nt  ly  a  v*r«Jon  of  the  hybrid  method 
has  be*n  sp*eiriraj|y  d*velopnd  at  DAe  to  estimate  the 
vaf1sii.»n  of  9*ro-)irt  pitching  moment  with  Mach  number. 

This  technlqu*  has  h*''n  un*d  l«i  a  design  process  to 

trim  drag  of  a  combat  aircraft  configuration, 
yielding  a  r-d»irtton  In  H  ft-dependent  drag  at  the 
crltlral  subno»i|r  sod  supersonic  d*slgn  points. 


|y i?C  fof  forebodlcs  and  pltot-lntake 

C  ow  i  s 

Th*  nA*  Fnlef  eofip  f<5f.  axlsymmetrlc  forebodlea^^ 
ban  usAd  »o  mlruiat*  t»»»  variation  of  drag  with 

Msai,  number  nf  Inw  mipersonlc  Speeds  for  the  forebody  of 
Fig  7?  at  Incidence.  Fig  ??  shows  that  the  method 

provide*  s  reas-^finbly  falthfuJ  representat Ion  of  the 
varla^lnn  f^r  MsaP  nuoihAns  h*tw**n  ]  and  1.?. 


iupmrMfile  Drag  Pwimiiit  Diit  to 
Wing  Coiwbor  M  TwttI 

1  mapfWB  —Bn  wp. 


As  n-'ted  In  seetion  A,  calculations  have  been  made  of  cow)  pressure  drag  by  a  ver-* 
slnti  nf  \hp  pA*  Fuier  cude?A  for  pitot  cowls  (Pig  ?••)  at  supersonic  speeds.  Fig  25  shows 
that  predtrfinns  by  title  m*thM  are  In  good  agreement  with  measurement. 

^  Eu  J  er/MuH  1  hJ  pek  joeUjod 

Althotigh  methods  auk'll  as  those  described  in  section  ft.l  and  A.?  have  demonstrated 
their  usefulttess  as  englne«rlng  tools,  Increasing  use  will  be  made  In  the  future  of 
methods  finch  as  th*  AHA/PAe  Eu 1 •r/Mul 1 1  block  code,  as  noted  in  section  The  appli¬ 

cation  or  thin  m*ih''d  to  wlng/body  configurations  representative  of  auperaonlr  combat 
aircraft  In  d*nrr1b*d  nnd  ass*ss*d  In  n*r  56.  In  this  study,  drag  ia  determined  by  the 
•joral’  and  ihun  n*«ds  to  b*  regarded  with  caution  becauae  of  the  sensitivity  of 

th*  method  to  disAfAtisstton  errors.  A  study  has  been  made  of  the  effect  on  drag  of  grid 
atnictui*  and  d*nnlty  but  this  wss  not  concluaiveis.  Therefore  the  aaaeaament  of  the 
w*th*d  ban  h*Aii  bssAd  mainly  upon  comparisons  with  measurements  of  wing  pressure  distri¬ 
bution*  and  overall  forAes  ms*l*  on  two  half  siodels.  In  order  that  the  comparison  Is  not 
srr*ct*d  hy  *Ki  ••aiiAoija  Af  frets,  such  ss  thes*  due  to  the  Interaction  between  the  hair  body 
and  tti*  sidewall  h'uindsry  Jay*r,  overall  force  messurementa  on  the  body  alone  are 
aubtrseted  (r'Mn  th*r.*  of  th*  wlnf'body  configuration  at  each  angle  of  incidence  and  an 
analogous  proce/lur*  la  used  in 
the  calculation  by  the  CPP 
method.  r.owparl ao»ts  are  shown 
In  Fig  3n  for  H  •  1.6  and  for 
one  of  the  wings  studied,  the 
calculated  value  of  drag  coef¬ 
ficient  having  been  Increased  by 
O.O05A  to  allow  for  akin  friction 
(aaaumed  to  be  unaffected  by  wing 
incidence,  thickness  and  camber). 

The  agreement  between  calculation 
and  sieaaurement  Is,  on  the  who!*, 
fair.  Pirferences  between  pre¬ 
diction  and  s*easiir*meMt  of  lift  at 
angles  of  incidence  above  about  6* 
can  be  *Kplaln*d  by  th*  effects  of 
shock-induced  separation  not  almu- 
lat*d  In  th*  cnlculatlori  method. 

The  obvloua  discrepancies  between 
calculation  and  measurement  of  drag 

at  low  lift  Is  believed  to  he  doe  mainly  to  Inaccurate  predictions  of  auction  near  the 
leadl»ig-*dge. 


flf  It  0««rflt  Sfl.  Sr»|  wil  sStMuy  waarnl  tterwi-nn  t»twr»s  ISvwy 
•et  ftartmeBl.  h  »■  1*.  Wls|  A 
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•  •  5  lUll  *s  multi g rid  tinethod 


Woodward'’  hRj^  uasd  Mall'e  niultlgrld  laethod^^  for  tolvlng  the  Euler  equation,  pre- 
viouRly  mentlon<*d  !n  eectlon  3.1.  to  etudy  the  wave  drag  of  aerofolla  with  rounded 
leadlng-»d*»a  at  aup^raonic  free  ntreaei  epeeda.  Thla  wethod  Is  particularly  suitable  for 
etuiylna  flows  of  tliJs  kind  alnce  It  has  an  unusually  large  nunber  of  grid  points  In  the 
loading-^dg^  region  and  la  thua  able  to  represent  accurately  the  strong  detached  shock  and 
t.h*  rapid  spatial  r|,«nKoa  In  the  flow  near  the  leading  edge. 


^Ig  31  n  I'lal  rates  some  of 
the  results  obtained  by  Woodward 
for  wave  drag  by  the  ’iocal’ 
isethod  and  shows  th^  effect  on  the 
variation  of  wave  drag  with  lift 
of  changing  nose  radius.  At  sero 
lift  an  optlmnni  nose  radius  of 
about  1|S  chord  is  obtained  hut, 
as  lift  Increases,  the  optimum 
value  becomes  smaller.  This 
Interesting  result.  Illustrates 
well  the  ability  of  CFO  to  provide 
re1ft^lv•ly  rapid  an.nrRjsments  of 
drag  dirferenreo  dti«  tn  changes  In 
ahap«  nnd  » he  means  or  determining 
drag  ep^  _ 


''M||nt,M|'in)5  njvMAnK.'^ 


’Vh  1  a  pnpei-  •^h-'vin  lhat, 

while  Ibr  wh*']lY  1.  Iie-.rel  1  ra  1  pre- 
d1'*i!e,,  of  aircraft  rlrag  to  net 

yet  persJble^  I'Fp  meljtofis  ejtot  In 
•  fK  r>'r  firsg  j  r'''1i''i  to. I  wt»1''li  ar® 
of  ent>«  t  de  ra  b  1  e  valiio  to 

alffiaft  designer  if,  t-|,c  rojlewlng 

1  act...; 


7»r«  Nil 


rtSkii 


)t  (•ItuUldf  «srl«lton  of  w«v«  drag  wtfh  lifi  for  »tr|rk>g 
nose  •  rsdhii  sorofoilt.  H  s  1.1 


.  5^10.. I  I.. II  -M  iiir.  fiiiap^.  or  .Ircrnf.  oompon.iit..  nt  th.  nr.ltmln.rv  etait»s  of  th» 

•,or.  Ipn; 

-  oii-iiy«,*i  of  H,ap  (»M-'  d  I  ap,ti'>. ,  a  oT  Bonre.  of  unw.nt.d  drag; 

-  ’oi I  rupo ( If  (oi, '  of  wfiid-l'inn.l  drag  data  to  TuJt  aca}a'. 

Fill'll, or  I  .1  1 1,0„,0|,|  ,  flp,  II.OI,.,,  .n  numarlpal  mathodd  for  solving  Hie  F.uler  eijuatlon 
to  rediif*  Mio  poiioifivity  or  drag  predlrtlona  Oy  these  methods  to  grid  density.  Such 
dfVflopsioiii's  would  silos  mnlUMfifk  arhamee  to  h*  esplolted  to  csleulate  the  drsg  or 
compiaj  eoiirig.iiaHoiis,  and,  as  such,  sould  be  a  step  In  the  direction  away  from  the 
current  dependence  n.,  wind-tunnel  tests. 

iiK  meiiiods  ef  solving  the  Reynolds-sversgeif ,  H«vler«Stckes  equation  have  yet  to 
mat  ®  a  signUicani  Impart  as  techniques  for  drsg  prediction.  Future  developments  in  this 
nr®a  depend  mainly  en  Improvements  being  mart*  to  the  turbulence  aiodels  used,  and  the 
pronpeei.fl  of  there  being  effecta'l  In  the  nesr  term  *re  uncertain.  Thus  vlscous/lnvtacXd 
lnier-sei in„  teehniqiier  npe  eipected  to  continue  to  feature  prominently  In  UK  drag  pre- 
diefien  methods  fcir  some  time  to  come. 
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